sum 总现值, detail
egen mode=mode(总现值)


************人均价值加房贷***************
replace 房贷=人均价值 if 总现值> 500000

replace 房贷=房贷价值 if 总现值<= 500000

egen sum_房贷=sum(房贷) 


egen sum_房贷1=sum(房贷) if 收入阶级==1
gen 税负1=sum_房贷1/sum_房贷 if 收入阶级==1

egen sum_房贷2=sum(房贷) if 收入阶级==2
gen 税负2=sum_房贷2/sum_房贷 if 收入阶级==2

egen sum_房贷3=sum(房贷) if 收入阶级==3
gen 税负3=sum_房贷3/sum_房贷 if 收入阶级==3

egen sum_房贷4=sum(房贷) if 收入阶级==4
gen 税负4=sum_房贷4/sum_房贷 if 收入阶级==4

***********公共服务均等分配*************
egen sum_number=sum(number)
gen 人均房贷福利=sum_房贷/sum_number
gen 家庭房贷福利=number*人均房贷福利
gen 房贷净福利=家庭房贷福利-房贷
sum 房贷净福利 if 房贷净福利>0
sum 房贷净福利 if 房贷净福利<=0
sum 总收入 if 房贷净福利>0
sum 总收入 if 房贷净福利<=0
sum 自有套数 if 房贷净福利>0
sum 自有套数 if 房贷净福利<=0
sum number if 房贷净福利>0
sum number if 房贷净福利<=0
sum 人均建筑面积 if 房贷净福利>0
sum 人均建筑面积 if 房贷净福利<=0
sum 总房贷 if 房贷净福利>0
sum 总房贷 if 房贷净福利<=0
egen sum_房贷净福利=sum(房贷净福利)
sum 房贷净福利 
***********仅分配给有未成年子女的家庭*************
egen sum_房贷=sum(房贷)
egen sum_未成年=sum(未成年)
gen 人均未成年房贷福利=sum_房贷/sum_未成年
gen 家庭未成年房贷福利=人均未成年房贷福利*未成年
gen 未成年房贷净福利=家庭未成年房贷福利-房贷
sum 未成年房贷净福利 if 未成年房贷净福利>0
sum 未成年房贷净福利 if 未成年房贷净福利<=0
sum 总收入 if 未成年房贷净福利>0
sum 总收入 if 未成年房贷净福利<=0
sum 自有套数 if 未成年房贷净福利>0
sum 自有套数 if 未成年房贷净福利<=0
sum number if 未成年房贷净福利>0
sum number if 未成年房贷净福利<=0
sum 人均建筑面积 if 未成年房贷净福利>0
sum 人均建筑面积 if 未成年房贷净福利<=0
sum 总房贷 if 未成年房贷净福利>0
sum 总房贷 if 未成年房贷净福利<=0
egen sum_未成年房贷净福利=sum(未成年房贷净福利)
sum 未成年房贷净福利

************净福利决定因素********************
gen 总面积=一套面积+二套面积+三套面积+四套面积+五套面积+六套面积
replace 是否拆迁=0 if 是否拆迁==2
replace 住房=1 if j1001==1
sum  房贷净福利 number 总收入 总现值 自有套数 人均建筑面积  总储蓄 住房 是否拆迁 房价 税率
gen ln房贷净福利=ln(房贷价值净福利)
gen ln总收入=ln(总收入)
gen ln总现值=ln(总现值)
gen ln总储蓄=ln(总储蓄)
gen ln房价=ln(房价)
gen ln人均建筑面积=ln(人均建筑面积)
gen 无减免福利=-房贷
sum 无减免福利
sum 房贷净福利
sum 未成年房贷净福利
sum 税率
reg 无减免福利 number 自有套数 总面积 总现值 总收入 总储蓄  基础消费  住房 是否拆迁 税率, r 
reg 房贷净福利 number 自有套数 总面积 总现值 总收入 总储蓄  基础消费  住房 是否拆迁 税率, r 
reg 未成年房贷净福利 number 自有套数 总面积 总现值 总收入 总储蓄 基础消费 住房 是否拆迁 税率, r

preserve
keep if 省份1==31 | 省份1==32 | 省份1==33 | 省份1==35 | 省份1==37 | 省份1==11 | 省份1==12 | 省份1==13 | 省份1==44 | 省份1==46 | 省份1==21 | 省份1==22 | 省份1==23
reg 未成年房贷净福利 number 自有套数 总面积 总现值 总收入 总储蓄 基础消费 住房 是否拆迁 税率, r
restore 

preserve
keep if 省份1==14 | 省份1==34 | 省份1==36 | 省份1==41 | 省份1==42 | 省份1==43 | 省份1==45  
reg 未成年房贷净福利 number 自有套数 总面积 总现值 总收入 总储蓄 基础消费 住房 是否拆迁 税率, r
restore 

preserve
keep if 省份1==61 | 省份1==62 | 省份1==63 | 省份1==64 | 省份1==50 | 省份1==51 | 省份1==52 | 省份1==53
reg 未成年房贷净福利 number 自有套数 总面积 总现值 总收入 总储蓄 基础消费 住房 是否拆迁 税率, r
restore 


tab 房价
preserve
keep if 房价>7000
reg 未成年房贷净福利 number 自有套数 总面积 总现值 总收入 总储蓄 基础消费 住房 是否拆迁 税率, r
restore 

preserve
keep if 房价<=7000 
reg 未成年房贷净福利 number 自有套数 总面积 总现值 总收入 总储蓄 基础消费 住房 是否拆迁 税率, r
restore 


replace 税后房贷指数1=总收入-房贷  //房贷价值
sum 税后房贷指数1
replace 税后房贷指数1=税后房贷指数1/104950.4 
replace 税后房贷指数1=税后房贷指数1^(0.5)  //系数为0.5 
egen sum_税后房贷指数1=sum(税后房贷指数1)
replace 税后房贷指数1=sum_税后房贷指数1/8727
replace 税后房贷指数1=税后房贷指数1^2
replace 税后房贷指数1=1-税后房贷指数1

replace 税后房贷指数2=总收入-房贷  //房贷价值
replace 税后房贷指数2=log(税后房贷指数2) //系数为1
egen sum_税后房贷指数2=sum(税后房贷指数2)
replace 税后房贷指数2=sum_税后房贷指数2/8727
replace 税后房贷指数2=2.71828^(税后房贷指数2)  //ye
replace 税后房贷指数2=1-(税后房贷指数2/104950.4)

replace 税后房贷指数3=总收入-房贷  //房贷价值
replace 税后房贷指数3=税后房贷指数3/104950.4 
replace 税后房贷指数3=税后房贷指数3^(-1)  //系数为2
egen sum_税后房贷指数3=sum(税后房贷指数3)
replace 税后房贷指数3=sum_税后房贷指数3/8727
replace 税后房贷指数3=税后房贷指数3^(-1)
replace 税后房贷指数3=1-税后房贷指数3

replace 税后房贷福利指数1=总收入+房贷净福利  // 家庭房贷价值福利
sum 税后房贷福利指数1
replace 税后房贷福利指数1=税后房贷福利指数1/106647.8
replace 税后房贷福利指数1=税后房贷福利指数1^(0.5)  //系数为0.5 
egen sum_税后房贷福利指数1=sum(税后房贷福利指数1)
replace 税后房贷福利指数1=sum_税后房贷福利指数1/8727
replace 税后房贷福利指数1=税后房贷福利指数1^2
replace 税后房贷福利指数1=1-税后房贷福利指数1

replace 税后房贷福利指数2=总收入+房贷净福利  //家庭房贷价值福利
replace 税后房贷福利指数2=log(税后房贷福利指数2) //系数为1
egen sum_税后房贷福利指数2=sum(税后房贷福利指数2)
replace 税后房贷福利指数2=sum_税后房贷福利指数2/8727
replace 税后房贷福利指数2=2.71828^(税后房贷福利指数2)  //ye
replace 税后房贷福利指数2=1-(税后房贷福利指数2/106647.8)


replace 税后房贷福利指数3=总收入+房贷净福利  // 家庭房贷价值福利
replace 税后房贷福利指数3=税后房贷福利指数3/106647.8
replace 税后房贷福利指数3=税后房贷福利指数3^(-1)  //系数为2
egen sum_税后房贷福利指数3=sum(税后房贷福利指数3)
replace 税后房贷福利指数3=sum_税后房贷福利指数3/8727
replace 税后房贷福利指数3=税后房贷福利指数3^(-1)
replace 税后房贷福利指数3=1-税后房贷福利指数3

replace 税后未成年福利指数1=总收入+未成年房贷净福利  // 家庭未成年房贷价值福利
sum 税后未成年福利指数1
replace 税后未成年福利指数1=税后未成年福利指数1/106647.8
replace 税后未成年福利指数1=税后未成年福利指数1^(0.5)  //系数为0.5 
egen sum_税后未成年福利指数1=sum(税后未成年福利指数1)
replace 税后未成年福利指数1=sum_税后未成年福利指数1/8727
replace 税后未成年福利指数1=税后未成年福利指数1^2
replace 税后未成年福利指数1=1-税后未成年福利指数1

replace 税后未成年福利指数2=总收入+未成年房贷净福利  //家庭未成年房贷价值福利
replace 税后未成年福利指数2=log(税后未成年福利指数2) //系数为1
egen sum_税后未成年福利指数2=sum(税后未成年福利指数2)
replace 税后未成年福利指数2=sum_税后未成年福利指数2/8727
replace 税后未成年福利指数2=2.71828^(税后未成年福利指数2)  //ye
replace 税后未成年福利指数2=1-(税后未成年福利指数2/106647.8)

replace 税后未成年福利指数3=总收入+未成年房贷净福利  // 家庭未成年房贷价值福利
replace 税后未成年福利指数3=税后未成年福利指数3/106647.8
replace 税后未成年福利指数3=税后未成年福利指数3^(-1)  //系数为2
egen sum_税后未成年福利指数3=sum(税后未成年福利指数3)
replace 税后未成年福利指数3=sum_税后未成年福利指数3/8727
replace 税后未成年福利指数3=税后未成年福利指数3^(-1)
replace 税后未成年福利指数3=1-税后未成年福利指数3

replace 社会福利房贷福利指数1=((106647.8 *(1-0.375))/(106647.8*(1-0.386)))^(-1)  //系数为0.5
gen 社会福利房贷福利指数2=(109897.3 *(1-0.0049))/(109101.6*(1-0.0044))  //系数为1
gen 社会福利房贷福利指数3=((109897.3  *(1-0.0097))/(109101.6*(1-0.0088)))^(-1)  //系数为2

sum 税后未成年福利指数3 if 收入阶级==4 

sum 总收入 if 收入阶级==4
gen 税前指数1=总收入/109101.6   if 收入阶级==2
replace 税前指数1=税前指数1^(0.5) if 收入阶级==2 //系数为0.5 
egen sum_税前指数1=sum(税前指数1) if 收入阶级==2
replace 税前指数1=sum_税前指数1/ 2133  if 收入阶级==2
replace 税前指数1=税前指数1^2 if 收入阶级==2
replace 税前指数1=1-税前指数1 if 收入阶级==2

gen 税前指数2=log(总收入)  if 收入阶级==4 //系数为1
egen sum_税前指数2=sum(税前指数2) if 收入阶级==4
replace 税前指数2=sum_税前指数2/669  if 收入阶级==4
replace 税前指数2=2.71828^(税前指数2) if 收入阶级==4 //ye
replace 税前指数2=1-(税前指数2/ 347802.4 ) if 收入阶级==4

gen 税前指数3=总收入/  109101.6    if 收入阶级==2
replace 税前指数3=税前指数3^(-1) if 收入阶级==2 //系数为2
egen sum_税前指数3=sum(税前指数3) if 收入阶级==2
replace 税前指数3=sum_税前指数3/ 2133   if 收入阶级==2
replace 税前指数3=税前指数3^(-1) if 收入阶级==2
replace 税前指数3=1-税前指数3 if 收入阶级==2


gen 税后房贷福利指数1=总收入+未成年房贷净福利  if 收入阶级==4 // 家庭未成年房贷价值福利
sum 税后房贷福利指数1 if 收入阶级==4 
replace 税后房贷福利指数1=税后房贷福利指数1/337642    if 收入阶级==4
replace 税后房贷福利指数1=税后房贷福利指数1^(0.5) if 收入阶级==4 //系数为0.5 
egen sum_税后房贷福利指数1=sum(税后房贷福利指数1) if 收入阶级==4
replace 税后房贷福利指数1=sum_税后房贷福利指数1/669  if 收入阶级==4
replace 税后房贷福利指数1=税后房贷福利指数1^2 if 收入阶级==4
replace 税后房贷福利指数1=1-税后房贷福利指数1 if 收入阶级==4

gen 税后房贷福利指数2=总收入+未成年房贷净福利 if 收入阶级==4 //家庭未成年房贷价值福利
replace 税后房贷福利指数2=log(税后房贷福利指数2) if 收入阶级==4 //系数为1
egen sum_税后房贷福利指数2=sum(税后房贷福利指数2) if 收入阶级==4
replace 税后房贷福利指数2=sum_税后房贷福利指数2/ 669  if 收入阶级==4
replace 税后房贷福利指数2=2.71828^(税后房贷福利指数2) if 收入阶级==4 //ye
replace 税后房贷福利指数2=1-(税后房贷福利指数2/337642 ) if 收入阶级==4

gen 税后房贷福利指数3=总收入+未成年房贷净福利 if 收入阶级==4 // 家庭未成年房贷价值福利
replace 税后房贷福利指数3=税后房贷福利指数3/337642     if 收入阶级==4
replace 税后房贷福利指数3=税后房贷福利指数3^(-1) if 收入阶级==4 //系数为2
egen sum_税后房贷福利指数3=sum(税后房贷福利指数3) if 收入阶级==4
replace 税后房贷福利指数3=sum_税后房贷福利指数3/ 669  if 收入阶级==4
replace 税后房贷福利指数3=税后房贷福利指数3^(-1) if 收入阶级==4
replace 税后房贷福利指数3=1-税后房贷福利指数3 if 收入阶级==4


gen 社会福利房贷福利指数1=((337642.0 *(1-0.0832))/(347802.4*(1-0.0860)))^0.5  if 收入阶级==4 //系数为0.5
gen 社会福利房贷福利指数2=(337642.0 *(1-0.1317))/(347802.4*(1-0.1366))  if 收入阶级==4 //系数为1
gen 社会福利房贷福利指数3=((337642.0 *(1-0.1843))/(347802.4*(1-0.1913)))^(-1)  if 收入阶级==4 //系数为2